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ABSTRACT 

We carry out a comprehensive analysis of the behavior of the magnetorotational instabihty (MRI) in viscous, 
resistive plasmas. We find exact, non-linear solutions of the non-ideal magnetohydrodynamic (MHD) equations 
describing the local dynamics of an incompressible, differentially rotating background threaded by a vertical 
magnetic field when disturbances with wavenumbers perpendicular to the shear are considered. We provide a 
geometrical description of these viscous, resistive MRI modes and show how their physical structure is modified 
as a function of the Reynolds and magnetic Reynolds numbers. We demonstrate that when finite dissipative 
effects are considered, velocity and magnetic field disturbances are no longer orthogonal (as it is the case in 
the ideal MHD limit) unless the magnetic Prandtl number is unity. We generalize previous results found in the 
ideal hmit and show that a series of key properties of the mean Reynolds and Maxwell stresses also hold for 
the viscous, resistive MRI. In particular, we show that the Reynolds stress is always positive and the Maxwell 
stress is always negative. Therefore, even in the presence of viscosity and resistivity, the total mean angular 
momentum transport is always directed outwards. We also find that, for any combination of the Reynolds and 
magnetic Reynolds numbers, magnetic disturbances dominate both the energetics and the transport of angular 
momentum and that the total mean energy density is an upper bound for the total mean stress responsible for 
angular momentum transport. The ratios between the Maxwell and Reynolds stresses and between magnetic 
and kinetic energy densities increase with decreasing Reynolds numbers for any magnetic Reynolds number; 
the lowest limit of both ratios is reached in the ideal MHD regime. 

Subject headings: black hole physics — accretion, accretion disks — MHD — instability — turbulence 



1. INTRODUCTION 

The magnetorotational instability (MRI, iBalbus & HawlevI 
1199 iL 11998 ) has been widely studied in the inviscid and 
perfectly conducting, magnetohydrodynamic (MHD) limit. 
The departures from this idealized situation are usually 
parametrized according to the Reynolds Yie — vl/v and mag- 
netic Reynolds Rm = vl/rj numbers, where v and I stand 
for the relevant characteristic velocity and lengthscale and v 
and rj stand for the kinematic viscosity and resistivity. The 
ideal MHD regime is then formally identified with the limit 
Re , Rm oo. There are many situations of interest in 
which the effects of dissipation need to be considered. 

From the astrophysical point of view, accretion disks 
around young stellar objects constitute one of the most com- 
pelling reasons for investigating the MRI beyond the ideal 
limit. In particular, there is great interest in understanding 
to what extent can MHD turbulence driven by the MRI en- 
able efficient angular momentum transport in cool, poorly 
conducting , protoplanetary disks (s ee, e.g. iBlaes & BalbusI 
1991 ljiiT[l99l iGammid [T996I: [SanoA Mivama 1992 
Salmeron & Wardld 120051) Most of the studies addressing 
the effects of dissipation in non-ideal MRI have usually fo- 
cused in inviscid, resistive plasmas. However, accretion disks 
are characterized by a wide range of magnetic Prandtl num- 
bers, with Pm ^ V jj] varying by several orders of magnitude 
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across the entire disk (see, e.g.. 'Balbu s & Henrill2007 l). In 
order to understand the behavior of the MRI under these con- 
ditions it is necessary to relax the assumption of an inviscid 
plasma. 

A large fraction of the shearing box simulations address- 
ing the non-linear regime of the MRI have been carried out in 
the ideal MHD limit, i.e., without including explicit dissipa- 
tion in the codes (see, e.g..[Hawlev. Gammie. & Balbus 1995} 
[Brandenburg et ak 1995; Sano et al. 2004). However, even in 
the absence of explicit viscosity and resistivity, finite differ- 
ence discretization leads to numerical diffusion/dissipation. 
Therefore, even in this type of simulations, it is necessary to 
understand the impact of these numerical artifacts that lead to 
departures from the ideal MHD regime and how similar they 
are when compared with physical (resolved) dissipation. 

A handful of numerical studies with explicit resistivity but 
zero physical viscosity have been carried out in order to un- 
derstand the effects of ohmic dissipation in the saturation of 
MRI-d riven turbulence ( s ee, e.g .,p.ano. Inutsuka, & Mivama 
1998"; 'Sano & Inutsuka' 1200 U iFl eming, Stone, & Hawlev 
2000; Sano et al. 2004; T urner. Sano . & Dziourkevitch2007). 
In particular, iSano & Ston3 (12003) have shown that the sat- 
uration level of the stresses increases with increasing mag- 
netic Reynolds number and seem to converge to an asymptotic 
value for magnetic Reynolds numbers larger than unity. 

Recent work has pointed out problems with conver- 
gence in zero-net-flux numerical simulations of ideal 
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MHD driven by the MRI (iPessah. Chan. & Psalfel 120071: 
iFromang & Papaloizoul200'7l) implying the necessity of incor- 
porating expHcit dissipation in the codes. Numerical studies 
with both resistivity and v iscosity, in the p r esenc e of a mean 
vertical magnetic field (Lesur & Longaretti 2007) and in the 
case of zero net flux (Fromang et al- 2007), have begun to un- 
cover how the characteristics of fully developed MRI-driven 
turbulence depends on the Reynolds and magnetic Reynolds 
numbers. Even though the ranges in Reynolds and magnetic 
Reynolds numbers that can be currently addressed is still lim- 
ited, the results obtained from the simulations suggest that the 
magnetic and kinetic energies contained in turbulent motions 
in the saturated regime depend on the values of the microphys- 
ical viscosity and resistivity. In particular, the mean angular 
momentum transport in the turbulent state increases with in- 
creasing magnetic Prandtl number. 

From the experimental perspective, understanding the ef- 
fects of non-vanishing resistivity and viscosity in the behav- 
ior of the MRI seems imperative, since the physical condi- 
tions achievable in the la boratory depart significantly from 
the ideal MHD regi me (|Ji. G oodman & Kagevama' '2001!: 
Goodman & Ji 200 2; ISisaneTa l. 2004; Liu, Goodman, & Ji 
2006c iRudiger. Schultz. Shalvbkovl 12003ft . Liquid metals 
(such as sodium, gallium, and mercury) are often character- 
ized by rather low magnetic Prandtl numbers (Pm ~ 10^^- 
IQ-'^). Although the regime of Reynolds numbers involved is 
still orders of magnitude smaller than any astrophysical sys- 
tem with similar magnetic Prandtl numbers, MRI experiments 
offer one of the few prospects of studying anything close to 
MHD astrophysical processes in the laboratory. 

A number of analyses addressing some aspects of the im- 
pact of viscosity and resistivity on the MRI in various dis- 
sipative limits appear scattered throughout the literature on 
theoretical, numeri cal, and experimental MRI. More recently, 
iLesaffre & BalbusI ([2007) have found particular solutions of 
the viscous, resistive MHD equations (including even a cool- 
ing term) in the shearing box approximation. However, we 
are unaware of any comprehensive, systematic study address- 
ing how the MRI behaves in viscous, resistive, differentially 
rotating magnetized plasmas for arbitrary combinations of the 
Reynolds and magnetic Reynolds numbers. The aim of this 
work is to carry out this analysis in detail. 

The rest of paper is organized as follows. In § |2l we 
state our assumptions. In §[3] we solve the eigenvalue prob- 
lem defined by the MRI for arbitrary Reynolds and mag- 
netic Reynolds numbers. We provide closed analytical ex- 
pressions for the eigenfrequencies and the associated eigen- 
vectors. In §2] we address the unexplored physical structure 
of MRI modes for finite Reynolds and magnetic Reynolds 
numbers and derive simple analytical expressions that de- 
scribe these modes in various asymptotic regimes. In § |6l 
we calculate the correlations between magnetic and veloc- 
ity MRI-driven perturbations that are related to angular mo- 
mentum transport and energy densities. We find that some 
key results previously shown to hold in the ideal MHD limit 
jPessah. Chan. & Psaltisl l2006) are also valid in the non-ideal 
regime. In particular, we show that even though the effective- 
ness with which the MRI disrupts the laminar flow depends on 
the Reynolds and magnetic Reynolds numbers, the instabil- 
ity always transports angular momentum outwards. We also 
find that magnetic perturbations dominate both the energetics 
and the transport of angular momentum for any combination 
of the Reynolds and magnetic Reynolds numbers. In §|7]we 



summarize our findings and discuss the implications of our 
study. 

2. ASSUMPTIONS 

Let us consider a cylindrical, incompressible background 
characterized by an angular velocity profile = rt{T)z, 
threaded by a vertical magnetic field B ~ B^z. We work in 
the shearing box approximation, which consist of a first order 
expansion in the variable r — ro of all the quantities charac- 
terizing the flow at the fiducial radius tq. The goal of this 
expansion is to retain the most relevant terms governing the 
dynamics of the MHD fluid in a locally-Cartesian coordinate 
system co-orbiting and corrotating with the background flow 
with local (Eulerian) velocity v — rg ^o<t). (For a more de- 
tailed discussion on this expansion see iGoodman & 5CulfT994l 
and references therein.) 

The equations governing the dynamics of an incompress- 
ible MHD fluid with constant kinematic viscosity ly and resis- 
tivity rj in the shearing box limit are given by 



dv 
'dt 



dB 



_ _V P + — + ^ ^ ' + vV^v , 

(1) 
(2) 



^ + {vV) B = (B- V) V + rjV'^B , 
where P is the pressure, p is the (constant) density, the factor 



q ■ 



d\ni 



(3) 



parametrizes the magnitude of the local shear, and we have 
defined the (locaUy-Cartesian) differential operator 



„ d <j} d d 
V = r h— V z — 

dr ro dcj) dz 



(4) 



where f, (f), and z are, coordinate-independent, orthonormal 
vectors corrotating with the background flow at tq. The con- 
tinuity equation reduces to V-t> — and there is no need for 
an equation of state since the pressure can be determined from 
this condition. 

We focus our attention on the dynamics of perturbations 
that depend only on the vertical coordinate. Under the current 
set of assumptions, these types of perturbations are known 
to exhibit the fastest growth rates in the ideal MHD case 
(Balbus & Hawley 1992, 1998; Pessah & Psaltis 2005). The 
equations governing the dynamics of these perturbations can 
be obtained by noting that the velocity and magnetic fields 
given by 

v = Svr{z)r + [~qno{r - ro) + Sv^(z)]cl) + 6vz{z)z ,(5) 
B = 5Br{z)f + 5B4z)^ + [B, + SB,{z)]z , (6) 

where the time dependence is implicit, constitute a family 
of exact, non-linear, solutions to the viscous, resistive MHD 
equations ([TJ-lO. As noted in Goodman & Xu ( 1994), even 
in the dissipative case, the only non-linear terms, which are 
present through the perturbed magnetic energy density, are ir- 
relevant in the case under consideration. 
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We can further simplify equations ^ and ^ by removing 
the background shear flow' f shear = — 9f^o(^ — fo)4> ^nd by 
realizing that we can take Svz{z) — 5Bz{z) — without loss 
of generality. We then obtain 



d_ 

dt 



Bz d 
A'Kp dz 



6Br 



— (5w0 = -(2 - q)noSv, 



92 
/ — 

dz' 



^ 47rp dz^ 



■Svr ■ 



(7) 



B. 



92 



dt 



5Br 



■Bz^SVr 

oz 



az^ 



SBr, 



(9) 



d 

—5B^ = -qna5Br 



B. 



92 



(10) 



'■dz ^^ ' ''dz^ 

where the first term on the right hand side of equation ([8]l is 
related to the epicyclic frequency 

K= ^2{2~q)no, (11) 
at which the flow variables oscillate in a perturbed hydrody- 
namic disk. For Keplerian rotation the parameter is g = 3/2 
and thus the epicyclic frequency is k = fio- 

It is convenient to define the new variables 6bi = 
SBi/^/4TTp for i — r,(t), and introduce dimensionless quan- 
tities by considering the characteristic time- and length-scales 
set by l/57o and Bz/{\/4:7rp flo). The equations satisfied by 
the dimensionless perturbations, 6vi, Sbi, are then given by 
dj:6Vr — 26v^ + dzSbr + ud^Svr , 
9j(5w0 = — (2 - q)5vr + diSb^ + I'djSv^ , 
ijdj6br , 
-qSbr + diSv^ + rjOf Sb^ , 

where t and z denote the dimensionless time and vertical co- 
ordinate, respectively. 

The dynamics of ideal MRI modes, with v = -q = 
0, is completely deterrn ined by the dimensionless shear q 
jPessah. Chan. & Psaltisll2006l) . The effects of viscous and 
resistive dissipation introduce two new dimensionless quan- 
tities that alter the characteristics and evolution of the MRI. 
With our choice of characteristic scales, it is natural to define 
the Reynolds and magnetic Reynolds numbers characterizing 
the MHD flow as ^ 

^ = 4, (16) 



dfSbr —dzSVr 



(12) 
(13) 
(14) 
(15) 



ReE 



RmE 



1 



(17) 



with associated magnetic Prandtl number 

„ Rm i> J/ 

Pm = = - = - . (18) 

Re fj rj 

In order to simplify the notation, we drop hereafter the tilde 
denoting the dimensionless quantities. In the rest of the paper, 
all the variables are to be regarded as dimensionless, unless 
otherwise specified. 

' In the shearing box approximation, the dependence of the background 
flow on the radial coordinate is strictly linear and therefore viscous dissipa- 
tion does not affect its dynamics. 

^ In a compressible fluid, the sound speed, Cg, provides another natural 
characteristic speed to define the Reynolds and magnetic Reynolds numbers. 
These defi nitio ns, e .g.. Re' and Rm', are related to those provided in equa- 
tions U6) and U7t via the plasma beta parameter /3 = (2/r)(cs/i'Az)^, 
with = Bz/\/4vrp the Alfven speed in the 2 direction, simply by 

Re' = (r/3/2)Re and Rm' = {r/3/2)Rm for a polytropic equation of 
state P = Kp^, with K and F constants. 



3. THE EIGENVALUE PROBLEM FOR THE NON-IDEAL MRI : 
A FORMAL ANALYTICAL SOLUTION 

In this section we provide a complete analytical solution 
to the set of equations (fT2]i-(fT5Tl as a function of the shear 
parameter q, (or, equivalently, the epicyclic frequency, k) for 
any set of values {ly, rj) defining the viscosity and resistivity. 

It is convenient to work in Fourier space, as this provides 
the advantage of obtaining explicitly the basis of modes that is 
needed to construct the most general solution satisfying equa- 
tions (fTSI i-lfTsll. Taking the Fourier transform of this set with 
respect to the ^-coordinate, we obtain the matrix equation 



dtS{kn,t) = LS{kn,t) 



where the vector d{kn,t) stands for 

8{kn,t) : 



SVr{kn,t) 

5vj{kn,t) 
5br{kn,t) 
Sb^{kn,t) 



(19) 



(20) 



and L represents the matrix 





-vkl 


2 




" 


L = 


-(2-9) 


-^kl 





ikji 









-Vkl 










ikji 


-q 


-r]kl . 



(21) 



The functions denoted by /(fc„, t) correspond to the Fourier 
transform of the real functions, /(z, <), and are defined via 



f{kn,t) = 



1 

2H 



H 



f{z,t)e-'''-'dz, 



(22) 



H 



where we have assumed periodic boundary conditions at z = 
±H, with H being the (dimensionless) scale-height and fc„ 
the wavenumber in the z-coordinate. 



_ nn 



(23) 



where n is an integer number. In order to simplify the nota- 
tion, hereafter we denote these wavenumbers simply by k. 

In order to solve the matrix equation ( fT9] l, it is convenient to 
find the eigenvector basis, {sj} with j — 1, 2, 3, 4, in which L 
is diagonal. This basis exists for all values of the wavenumber 
k (i.e., the rank of the matrix L is equal to 4, the dimension 
of the complex space) with the possible exception of a finite 
number of values of k. In this basis, the action of L over the 
set {ej} is equivalent to a scalar multiplication, i.e., 

idiag ej = (Tj for j = 1, 2, 3, 4 , (24) 

where {(Jj} are complex scalars. 

3.1. Eigenvalues 

In the eigenvector basis, the matrix L has a diagonal repre- 
sentation Ldiag = diag((Ti, (72, era, (74). The eigenvalues {(Tj}, 
with j ~ 1,2,3,4, are the roots of the characteristic poly- 
nomial associated with L, i.e., the dispersion relation associ- 
ated with the non-ideal MRI, which can be written in compact 
form as 



(fc2 + a^arjf + n^(k^ + al) - 4fc2 ^ , 



(25) 
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Fig. 1. — Growth rates 7+, eq. )44K as a function of the vertical wavenumber k for different combinations of Reynolds and magnetic Reynolds numbers for 
Keplerian rotation. In all three panels, the thick solid line corresponds to the ideal MHD limit, i.e., Ro, Rm — > 00. For any combination of the Reynolds and 
magnetic Reynolds numbers, the growth rate has a well defined, single maximum 7m ax that corresponds t o the most unstable mode fcmax- The range of unstable 
modes, < A; < fcc, is always finite, the critical wavenumber kc satisfies eq. (25} when cr = 0, see §|4J] Left: Growth rate 7+ for different values of the 
magnetic Reynolds number in the inviscid limit, i.e., Rc 00. The thin solid lines, in decreasing order, con'espond to Rm = 10, 1, 0.1. Middle: Growth 
rate 7+ for magnetic Prandtl number Pm = Rm/Rc = 1. The thin solid fines, in decreasing order, correspond to Re = Rm = 10, 1, 0.1. In all of the 
cases show n in the left and middle panels, the critical wavenumber, kc, below which unstable modes can exist decreases with increasing resistivity, see ij 14.31 
and i) |4.4| for analytic expressions of these marginally stable modes. Right: Growth rate 7+ for different values of the Reynolds number in the ideal conductor 
limit, i.e., Rm — > 00. The various curves, in decreasing order, correspond to Rc = 10, 1, . . . , 10 ~'^. I n this case, the range of unstable modes is insensitive to 
the Reynolds number, all the modes with wavenumbers shorter than fcc = v^2g are unstable, see § 14.21 It is evident that the growth rates and the characteristic 
scales, both fcmax and kc, are more sensitive to changes in the resistivity than to changes in the viscosity. The simultaneous analysis of all three panels leads 
to the conclusion that visco us, resistive modes with magnetic Prandtl number equal to unity resemble more closely inviscid, resistive modes rather than viscous, 
conductive ones, see ii l4.4l for the explanation of this behavior 



where we have defined the quantities 



and y is any of the solutions to the cubic equation 



0-„ = cr 



(26) 

j,^ — u -r rjk^ . (27) 

The dispersion relation ( |25l l is a fourth order polynomial with 
non-zero coefficients in a and . In order to find its roots 
it is convenient to take this polynomial to its depressed form. 
This can be achieved by defining the new variables and /i 
such that-' 



which has closed analytic solutions 

5 1 P 



aX 
"2" 



f3' 



with 



The resulting polynomial can then be written as 

/3cr„ + A = , 



acr,, 



where the coefficients a, (3, and A are given by 

a = 2{k^ - (/) + , 

A=(fc2_^2)2^^2(^2^^2- 

The solutions to equation (|30] | are 



4fc2 



4VA 



(28) 
(29) 

(30) 

(31) 
(32) 
(33) 

(34) 



0, (37) 

(38) 

(39) 
(40) 

(41) 



Note that the choice of either sign in U is immaterial. 

It is now trivial to write the solutions of the dispersion re- 
lation (125] ), cr, in terms of the variable a^. Using equations 
I we obtain 




cr = cr^ - i(l^ + ?7)P, 



(42) 



The eigenfrequencies of the viscous, resistive MRI modes are 
then given by 



where the subscripts a and b in the "+" and "— " signs label 
the four possible combinations of signs and we have defined 
the quantities'* 



2VA 



(43) 



3a y 

A = h -, 

4 2' 

A = {y + af-X, 



(35) 
(36) 



3.2. Four classes of solutions 



All of the quantities A, A, and y, depend on the viscosity 



and the resistivity rj only through /i~ 



{v - riY. This 



^ A physical interpretation of the variable ^ is provided in §[33] 
^ Defining the quantities A and A in this way allo ws u s to show explicitly 
that in the limit t', — > the solutions to equation {25\ converge smoothly 
to the solutions found in the ideal MHD case (see Appendix IaI. 



has a series of important implications, in particular, there is 
always a range of wavenumbers for which the discriminant in 
equation (|43] | is positive, i.e., VA — A > 0. It can also be seen 
that the last two terms between parentheses in equation ( l43T l 
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Fig. 2. — Left: Geometrical representation of the velocity field (black) and magnetic field (gray) perturbations for viscous, resistive MRI modes. Note that 
this is a projection of a single mode, which is inherently three-dimen sion al, onto the disk mid-plane (r, <t>,z = 0) . Th e velocity and magnetic field components 
are always out of phase in the vertical direction z by 7r/4, see eq. ^63^ The angles 9^ and 9^. defined in eqs. (50) and jSlj , re spec tively, correspond to the 
physical angles defining the planes (perpendicular to the disk midplane) containing the MRI-diiven perturbations, see eqs. (64) and (65). The relative magnitude 
of velocity and magnetic field perturbations is determined by eq. (54). Right: Evolution of the geometrical representation of the fastest-growing, non-ideal MRI 
mode, with associated wavenumber /cmax, with magnetic Prandtl number equal to unity, as a function of the Reynolds/magnetic Reynolds number. When the 
Reynolds/magnetic Reynolds number varies according to Re = Rm : oo ^ 0, the angles evolve according to : 7r/4 — > and d\, : 37r/4 — + 7r/2 and the 
relative amplitude of the perturbations evolves according to bo/vQ : 5/3 — > oo. Note that the velocity and magnetic field perturbations are always orthogonal 
for Pm = 1, see eq. (68). 

are always positive, i.e., \fK > and thus they always 

produce damping. Because of this, we can classify the modes 
in four types: two (damped) growing and decaying "unstable" 
modes with eigenvalues 

7±=±(VA-A)1/2 



1 - 



- ^fcM 1 



2" V' 2VA/ 2 V 2VA 
and two (damped) "oscillatory" modes with eigenvalues 



(44) 



2VA 



(45) 



2VA J 2 
We arbitrarily label these eigenvalues as 

(Ti = 7+ , (72 = 7- , era = iLu^ , (74 = iuj^ . (46) 

Figure [T] shows the growth rate 7+ as a function of the 
vertical wavenumber k for different combinations of the 
Reynolds and magnetic Reynolds numbers for Keplerian ro- 
tation. These growth rates are more sensitive to changes in 
the resistivity than to changes in the viscosity. A qualitative 
understanding of this behavior can be obtained by realizing 
that viscosity tends to quench the instability, without altering 
the large scale magnetic field. Thus, as long as the resistivity 
is negligible, the range of unstable lenghtscales are the same 
in both ideal and viscous, perfectly conducting fluids. On the 
other hand, resistivity tends to destroy the magnetic field at 
small scales having a stronger impact on the stability of the 
perturbations at these scales. 

Mathematically, the asymmetric response of the growth rate 
to changes in the viscosity i/ or the resistivity 77 originates in 
the different functional form of the terms that contribute to 



produce damping, i.e., (1 — k /2v A) and (1 + k /2\/ A), in 
the exponential growth characterized by 7+ in equation (l44l) . 
If the oscillatory modes, uj± in equation (05]), are considered 
instead, the roles of the plus and minus signs in these terms are 
interchanged. From this analysis we can infer that the "oscil- 
latory" mode is affected (damped) more strongly by viscosity 
than by resistivity. 

The simultaneous analysis of the various panels in Figure[T] 
leads to the conclusion that viscous, resistive unstable modes 
with magnetic Prandtl number equal to unity resemble more 
closely inviscid, resistive modes rather than viscous, conduc- 
tive ones. In §|4]we provide analytical expressions to support 
this conclusion. 

3.3. Normalized Eigenvectors: Geometrical Representation 

The set of normalized eigenvectors, {6^^ }, associated with 
the eigenvalues ( |46] l are given by 



where 



e,(fc) 



for J = 1,2,3,4, 



ik 



(47) 



— «/c[2(T,„ + q{i^ — r])k'^]/{k'^ + cr^jO-^j) 



(48) 

(Ti^j = (jj + vk^, (Tjjj — <jj + rjk^, and the norms are given by 



e, = 



.1=1 



1/2 



(49) 



Here, is the Z-th component of the (unnormalized) eigen- 
vector associated with the eigenvalue aj . 
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Fig. 3. — Critical wavenumber fcc, see eq. j70t . corresponding to the marginally stable MRI mode for Keplerian rotation in different dissipative regimes. The 
horizontal lines at fee = VS represent the ideal MHD limit, eq. i76\ . Left: Critical wavenumber fee as a functio n of the magnetic Reynolds number for different 
values of the Reynolds number The thick solid line denotes the inviscid limit, i.e., Rc oo. Note that eq. )87t des crib es this curve exactly. The thin solid 
hnes, in decreasing order, correspond to Re = 1, 0.1, . . .. For small magnetic Reynolds number fee Rm, see eq. )87t . For finite Reynolds numbers, such 
that ReRm < 1 there is a transition between the regimes Rm "C 1 and Rm 2> 1 such that fee oc Km^^^, see eq. (98). Middle: Critical wavenumber fee 
as a function of the Reynolds number for different magnetic Prandtl numbers. Pm increases/decreases by an order of magnitude for each curve to the left/light 
of the thick solid line denoting the Pm = 1 case. The dashed lines fee oc Re Pm = Rm are calculated according to eq. )87K which gives the correct result 
even for Pm > 1, provided that the Reynolds number is sufficiently small. The dotted lines fee oc (ReRm)^/^ are calculated according to eq. )98t . Right: 
Ciitical wavenumber fee as a function of the Reynolds number for different values of the magnetic Reynolds number The thick soHd line coiTesponds to the ideal 



conduct or l imit, i.e., Rm — > oo. The thin solid lines, in decreasing order, correspond to Re = 10"", 10 



, 10" 



see eq. <98K while for Reynolds numbers larger than a few, the ciitical wavenumber is independent of Re regardless of the value of Rm, 

The set of four eigenvectors {6^ }, together with the set of 
complex scalars {(Jj} in equation ( l46b . constitute the full so- 
lution to the eigenvalue problem defined by the MRI for any 
combination of the Reynolds and magnetic Reynolds num- 
bers. 



For small Reynolds numbers fee oc Re^/'^, 



where, for the sake of simplicity, we have omitted the sub- 
script j on the left hand side. The expressions for the nor- 
malized eigenvectors {e^ }, for j — 1,2,3, 4, are then given 
by 



A geometrical representation of the eigenvectors ( |47] ) can 



be brought to light by defining the angles i 
ing to 

tan( 



vj 



and I 



2a. 



tan 9\^j — — — — 



2cr. 



accord- 



(50) 



(51) 



Vq cos Uvj 

vq sin dyj 
ibo cos 6'bj 
ibo sin 



(55) 



It is interesting to note that in this geometric representation 
the dispersion relation dZSl ) can be obtained from the trigono- 
metric identity 

, „;„2 , 



COS 



= 1 . 



It is important to remark that each of the four eigenvectors 
define, in principle, four sets of angles {6vj,9hj}, for j = 
1,2,3,4. We label the angl es as sociated with the different 
types of modes discussed in ij |3.2| according to 

Ovl = , ^'v2 = (^v^ , ^v3 = ^ J ^v4 = " ' (52) 

with similar definitions corresponding to 6'bj for j — 1,2,3,4. 
Note that these angles are defined in spectral space and de- 
pend, in general, on the wavenumber k, the epicyclic fre- 
quency, K, the viscosity and the resistivity rj. The angles as- 
sociated with the modes labeled by 7+ and 7_ are always real 
while the ones associated with the modes w+ and ui^ are in 
general complex. For the sake of brevity, in what follows we 
will refer to the set of angles describing unstable MRI modes 
{0^+,^;^+} simply as {d,,e^}. 

A normalized version of the MRI eigenvectors can be ob- 
tained by multiplying the set of vectors in equation ( l48T l by 
the amplitudes 



where the expressions for 



cos 9v 



sm I 




2(j(fc2 



(56) 



(57) 



(58) 



can be obtained from the definition of the angle d^j in equa- 
tion jSOl l. 

3.4. Temporal Evolution 

In physical space, the most general solution to the set of 
equations (fT2] l-(fT5]l. i.e.. 



A, 



(53) 



where we have defined 

2kvQ 



evolves in time according to 



dVr{z, t) 
Sbr{z, t) 

Sb^{z, t) 



1 



{v - rTi)[Aar,j + q{p - r])k'^]k'^ 



1/2 



k 



d{k,t) e' 



where 



2(/c2 



(54) 



4 



(59) 



(60) 



(61) 
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Fig. 4. — Wavenumber fcmax corresponding to the fastest growing non-ideal MRI modes for Keplerian rotation in different dissipative regimes. The dot- 
dashed horizontal lines at fcmax = -^15/16 represent the ideal MHD limit, eq. (85). Left: Fastest growing mode fcmax as a function of the magnetic Reynolds 
number for different values of the Reynolds number. The thick solid line denotes the inviscid limit, i.e.. Re oo. The thin solid lines, in decreasing order, 
correspond to Re = 1, 0.1, ... , 10~^. For m agn etic Reynolds numbers larger than unity, this wavenumber is independent of Rm regardless of the value of Re. 
The dashed line, calculated according to eq. j90t , provides the connect asymptotic limit fcmax oc Rm for small magnetic Reynolds numbers. Middle: Fastest 
growing mode fcmax as a function of the Reynolds number for different values of the magnetic Prandtl number. From left to righ t, the curves correspond to 
Pm = 10^, 10^, . . . , 1 (thick solid line), . . . , 10"''. The dashed lines fcmax Re Pm — Rm are calculated according to eq. (90}, which leads to the coiTect 
result even for Pm > 1 provided that the Reynolds number is sufficiently small. The dotted line fcmax oc Re^^^ results from eq. UOU . Right: Fastest growing 
mode fcmax as a function of the Reynolds number for different values of the magnetic Reynolds number. The thick solid line coiresponds to the ideal conductor 
limit, i.e., Rm — > oo. The thin solid lines, in decreasing order, con'espond to Rm = 1, 0.1, 0.01. For Reynolds numbers larger than unity, the growth rate is 
independent of Re regardless of the value of Rm. The dotted line, calculated according to eq. (TOT), provides the coirect asymptotic limit fc 

max Rei/2 for 

small Reynolds number. 



with {<Tj} and {e^ }, for j = 1,2,3,4, given by equations 
(|46] | and dSSl l. The initial conditions a{k, 0) are related to 
the initial spectrum of perturbations, ^(fc,0), via a(fc,0) — 
Q^^ 6{k,0). Here, is the matrix for the change of coor- 
dinates from the standard basis to the normalized eigenvector 
basis^ and can be obtained by calculating the inverse of the 
matrix 

Q = [e^^ e^^ 60-3 e„^] . (62) 

The temporal evolution of a single MRI-unstable mode in 
physical space can be obtained from a linear combination of 
ea-j (fc) and Gg.^ (— fc) as defined in equation ( |55] l. In particular, 

setting ai(fc,0) = a^(— fc,0) = —i/V2m equation (l60l) and 
substituting the result in equation i6l[ we obtain 



tions in physical space, see Figure |2] with 

Sv^{z,t) 



tan0. 



tan 0bj = 



6vr{z, t) 
Sb(j,{z,t) 



= const. 



~ const. . 



6br{z, t) 

Finally, defining the angle 0bvj such that 

dhvj = Oh] ~ (^vj + 2 ) ' 
+ fi'bvj) 



which implies that tan 9\,j tan 
the fact that 



tan(6ii + 62) 



tan Oi + tan O2 
1 — tan Oi tan 62 



5{z,t) 



\/2( 



vq cos6'v sin(fcz) 
vq sin^v sin(fcz) 
feo cos6'b cos(fcz) 
60 sin 0b cos(fcz) 



it is not difficult to show that 



(63) 



tan 0b- 



vj 



(64) 
(65) 

(66) 
1, and using 

(67) 
(68) 



These solutions are of particular importance for the linear 
late-time evolution of MRI modes. Note that any reasonable 
spectrum of initial perturbations of the type used in numerical 
simulations of shearing boxes will have a non-zero compo- 
nent along the unstable eigenvector ■ If the value of the 
magnetic field is such that the MRI can be excited for given 
values of the viscosity and resistivity then the exponentially 
growing perturbations in physical space will evolve towards a 
mode of the form (l63T l dominated by the lengthscale fc = fcmax 
for which the growth rate reaches its maximum value 7max- 

Note that if a perturbation in physical space is composed by 
a single mode of the type described in 13.21 no matter which 
class, then the angles defined in equations (ISOl l and (ISTT i are 
constant in time and are identical to the physical angles be- 
tween the planes containing magnetic and velocity perturba- 

^ The eigenvectors (55) are not in general orthogonal, i.e., e^^ ■ e^^, ^ 

for jr ^ j' ■ If desired, an orthogonal basis can b e constructed using the Gr am- 
Schmidt orthogonalization procedure (see, e.g.. lHoffman & Kunzell97Th . 



This means that /i 7^ provides a measure of how non- 
orthogonal velocity and magnetic field perturbations are. 

It is evident that when the magnetic Prandtl number 
approaches unity viscous, resistive, MRJ-driven magnetic 
and velocity perturbations tend to be orthogonal, i.e.. 



tan t&nO 



iiff 



: — 1, and therefore 

6hj — 0\ 



TT 

2 



for Pm = 1 , 



(69) 



for every wavenumber fc. This is illustrated in Figure |2] 
which shows the evolution of the angles 0b and 6^ corre- 
sponding to the most unstable MRI mode as a function of 
the Reynolds/magnetic Reynolds number when the magnetic 
Prandtl number is equal to unity. 

4. PHYSICAL STRUCTURE OF MRI MODES 

The evolution of the physical structure of a single growing 
MRI mode with wavenumber < fc < fcc is characterized by 
its growth rate 7+, the relative magnitude between the ampli- 
tudes of magnetic and velocity field perturbations, bo /vq, and 
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Fig. 5. — Maximu m gr owth rate 7max for Keplerian rotation in different dissipative regimes. Tlie dot-daslied liorizontal lines at 7max = 3/4 represent the 
ideal MHD limit, eq. j82t . Left: Maximum growth rate as a function of the magnetic Reynolds number for different values of the Reynolds number The thick 
solid line denotes the inviscid limit, i.e.. Re — > oo. The thin solid lines, in decreasing order, con'espond to Re = 1, 0.1, ... , 10~ . For magn etic Reynolds 
numbers larger than unity, the growth rate is independent of Rm regardless of the value of Re. The dashed line, calculated according to eq. )9U , provides the 
correct asymptotic limit 7max oc Rm for small magnetic Reynolds numbers. Middle: Maximum growth rate 7max as a function of the Reynolds number for 
different values of the magnetic Prandtl number. From left to right the curves correspond to Pm = 10'', 10^, . . . , 1 (thick solid line), . . . , lO^''. The dashed 
lines 7max oc Re Pm = Rm are calculated according to eq. 19 IK which leads to the correct result even for Pm > 1 provided that the Reynolds number is 
sufficiently small. The dotted line 7max cc Re^''^ results from eq. U02t . Right: Maximum growth rate 7max as a function of the Reynolds number for different 
values of the magnetic Reynolds number The thick solid line corresponds to the ideal conductor limit, i.e., Rm oo. The thin solid lines, in decreasing order, 
correspond to Rm = 1, 0.1, 0.01. For Reynolds numbers larger than unity, the growth rate is independent of Re regardless of the value of Rm. The dotted line, 
calculated according to eq. )102t , provides the correct asymptotic limit 7max oc Re^^^ for small Reynolds number. 

the two angles defining the planes containing them, and 9^. 
For any reasonable spectrum of initial perturbations the mode 
that exhibits the fastest exponential growth, 7niax, which we 
refer to as fcmax, will dominate the dynamics of the late time 
evolution of the viscous, resistive MRI. It is therefore of par- 
ticular interest to characterize the physical properties of this 
fastest growing mode in different dissipative regimes. 



4. 1 . Marginal and Fastest Growing MRI-modes 

Because the eigenvalue associated with the unstable grow- 
ing mode, 7+, is always real for any combination of the 
Reynolds and magnetic Reynolds numbers, it is possible to 
find the marginally stable mode such that 7+(/cc) = 0. 
Setting cr = in equation (IZST l. we obtain a polynomial in fee 



0, 



(70) 



valid for any value of the viscosity and resistivity. Note that k^ 
sets the minimum domain height for numerical simulations of 
viscous, resistive MRJ-driven turbulence. Figure |3] shows the 
solutions of equation (iTOl i in various dissipative regimes for 
Keplerian rotation. The analytic solutions of equation dTOl l 
are algebraically complicated but their asymptotic limits are 
rather simple. We find expressions for this critical wavenum- 
ber in several regimes of interest below. 

In the ideal MHD limit, it is straightforward to find sim- 
ple analytical expressions for the most unstable wavenumber, 
fcmax, and its associated growth rate, 7max- However, the an- 
alytical expressions that we derived for the eigenfrequencies 
in the non-ideal case, equation (l46T l. are not amenable to the 
usual extremization procedure. More precisely, it is very chal- 
lenging to find the values of fc,nax and 7max that satisfy 



dk 



= 0. 



(71) 



Figures]?] and |5] show the solutions of this equation in various 
dissipative regimes. 

Another possible path to find the values of the wavenumber 
fcmax, and the associated growth rate, is to use the fact that 
7inax satisfies simultaneously the dispersion relation (IZST i and 



a polynomial in 7max- The largest of the roots of this poly- 
nomial is the desired maximum growth rate. It is possible 
to find fcmax following a similar methodology, but eliminat- 
ing between the two polynomials 7max instead. However, for 
arbitrary values of the viscosity and resistivity, both proce- 
dures lead to a seventh degree polynomial whose roots must 
be found numerically, defeating altogether the attempt to find 
analytical expressions for fcmax and 7niax- 

Using as a guide the results shown in Figures ]4] and ]5] we 
follow an alternative procedure. The goal is to find simple 
analytical expressions to describe the asymptotic behavior of 
the most unstable mode, fcmax, and the maximum growth rate, 
7max, in different dissipative regimes. It is evident from Fig- 
ure JT] that /cniax < 1 and 7inax < 1 for all the non-ideal MRI 
modes. This information can be used to simplify the disper- 
sion relation and its derivative so as to decrease their order 
without loosing vital information. This makes it possible to 
obtain manageable, but accurate, expressions for fcmax and 
7max in different limiting regimes. 

Figure ]6] shows contour plots for the critical wavenumber, 
fee, the most unstable wavenumber, /cmax, and the maximum 
growth rate, 7max, as a function of the Reynolds and magnetic 
Reynolds numbers for Keplerian rotation. In all three panels, 
lighter gray areas correspond to larger values of k^, fc^ax, and 
7max, respectively. Note that in all the cases, the functional 
form of the contours naturally divides the plane (Re, Rm) in 
three distinctive regions that we denote according to I (ideal), 
R (resistive), and V (viscous). Note that when the most unsta- 
ble wavenumber, /cmax, and the maximum growth rate, 7inax, 
are considered, these regions can be associated with the re- 
gions where Re, Rm ^ 1, Re ^ Rm, and Re ^ Rm, re- 
spectively. The overlap between these regions is not as clear 
when the critical wavenumber fee is considered and some care 
is needed when deriving approximated expressions for it. 

4.2. Ideal MRI Modes 

Let us first demonstrate briefly how the formalism pre- 
sented in § ]3]reduces to previously known results in the ideal 
MHD limit. In the absence of dissipation, the eigenvalues 
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Fig. 6. — Contour plots for the critical wavenumber, fee, the most unstable wavenumber, fcmax, and the maximum growth rate, 7max, for Keplerian rotation. In 
all three panels, Ughter gray areas correspond to larger values of fcc, fcmax, and 7max, respectively. The solid lines highlight the contours for kc = 1, . . . , 10~J 
and fcmax =^10_^, . . . , 10 ~*^. The labels I (idea l), R (resistive), and V (viscous), denote the thi'ee regions of the (Re, Rm) plane where equations (76), (85), 
and (82); 187) , j90) and )9U ; and (98), UOIK and )102) are valid, respectively. The dashed lines dividing the three regions are obtained by equating neighboring 
approximations for fc'c, fcmax, and 7max. 



{co,j}, with j = 1, 2, 3, 4, are the r oots of the dispersion re- 
lation associated with the ideal MRI jBalbus & HawlevI 19911 
[T998h . 

+ alj)"^ + + alj) - 4fc2 = , (72) 

and are given by dPessah. Chan. & Psaltisll2006l) 

^o,j =±(-Ao±yA;^j , (73) 
where we have defined the quantities Aq and Aq such that 

Ao 



Y 



(74) 
(75) 



The critical wavenumber for the onset of the ideal MRI is 
obtained by setting o-q = in the dispersion relation ( |72] i. this 
leads to 

fcc = = V4 - k2 . (76) 

For all the modes with wavenumbers k < kc the difference 
\/Ao — Aq is positive and we can define the "growth rate" 70 
and the "oscillation frequency" ojq by 



70 ; 




(77) 
(78) 



both of which are real and positive (for all positive values of 
the parameter q). This shows that two of the solutions of equa- 
tion ( |72] | are real and the other two are imaginary. We can thus 
write the four eigenvalues in compact notation as 



Co,! = 70 : 



Co, 2 



• 



(79) 



In the ideal MHD limit, it is evident that the velocity and 
magnetic field perturbations are orthogonal for any mode, i.e., 
tan^?v tan^b = ^ 1, see equation (|68] |. and therefore 



TT 

2 



(80) 



The temporal evolution of a single MRJ-unstable mode in 
physical space reduces to 



S{z,t) 



V2f 



vq cos 9v sin(fcz) 
Vq sin 0v sin(fcz) 
bo sin 6v cos(fcz) 
-bo cos 6v cos(fcz) 



(81) 



These are esse ntially the (normalized) perturbations found in 
equation (4) in lGoodman & Xul (1 19941) ^. 

From the definition of the angle 0~ 
be seen that 



It can 



see equation ( 

The maximum growth rate can be obtained by noting that 
Jo — Q sin 6'v cos 6^ and therefore the maximum growth cor- 
responds to 

7max ^2^4' 
It then follows that, in the absence of dissipation, the planes 
containing the exponentially growing velocity and magnetic 
field perturbations are characterized by the angles 

(83) 



4 ' 



(84) 



regardless of the value of the shearing parameter/epicyclic fre- 
quency. 

Finally, noting that the wavenumber for which the maxi- 
mum growth rate is realized is 



1 



16 



(85) 



and using equation ( l54b for the ratio between the amplitudes 
of the magnetic and velocity fields we obtain 



vo 



(86) 



In § 16.21 we derive equations for the MRI-driven Reynolds 
and Maxwell stresses, as well as the kinetic and magnetic en- 
ergy densities associated with the perturbations. Equations 



* Note that the angle 7 in lGoodman & Xj <1994l) . in our notation defined 
by tan 7 = —5Br/SB^, is such that 7 = in the positive azimuthal axis 
and it takes increasingly positive values in the counter-clockwise direction. 
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Fig. 7. — Opening angle, Sdiff = Sb ~ ^v, between the planes containing the fastest exponentially growi ng m agn etic and velocity perturbations for Keplerian 
rotation in various dissipative regimes. In the ideal MHD limit the opening angle is BdiS = tt/S, see eqs. (83) and )84t . Left: opening angle Sdift as a function 
of the magnetic Reynolds number for different values of the Reynolds number. The thick solid line denotes the inviscid limit, i.e.. Re — > oo. The thin solid 
lines, in decreasing order according to 9diff f™ fixed Rm, correspond to Re = 10, 1, . . . , 10~^^. For magnetic Reynolds numbers much larger than unity, the 
opening angle is independent of Rm regardless of the value of Re. For sufficiently small/large Rm, ddiS — > tt/Z provided that Re S> Rin. For sufficiently 
small/large Rm, Sdiff — > 7r/2 — arctan(K/2) provided that Rm <g; Re. Note that this coiresponds to 9diff = 63°26' for a Keplerian disk. The only conditions 
under which Sdift exceeds tt/2 are such that Re > Rm ~ 1. Middle: opening angle 9diff as a fu nction of the Reynolds number for different values of the 
magnetic Prandtl number. The thick solid line at 9diff = ''^/'^ coiresponds to Pm = 1, see eq. i69\ and Fig.|2] The thin solid lines with peaks at 9difl > ''"/2 
correspond, from left to right, to Pm = 10~^ , 10~^, .... The thin solid lines, with Sdiff < 7r/2, in decreasing order according to 9diff f™ fixed Re, correspond 
to Pm = 10, 10^, . . .. Note that for small Re, Sdiff ""72 for Pm ~ 1, while Sdiff '■"/2 — arctan(K/2) for Pm ^ 1. Right: opening angle 9diff 
as a function of the Reynolds number for different values of the magnetic Reynolds number. The thick solid line coiresponds to the ideal conductor limit, i.e.. 



Rm — > oo. The thin solid lines, from right to left, con'espond to Rm : 
of Re regardless of the value of Rm. 



10,. 



,10 . For Reynolds numbers larger than unity, the opening angle is independent 



(fT2TT l. (fT22] |. (fT29] l. and (flJOl l. show why, in the ideal MHD 
limit, equations (l83T l, (|84] |. and (|86] | are the reason for which 
the ratio between the Maxwell to the Reynolds stresses is 
identical to the ratio between magnetic and kinetic energy 
densities for any shear parameter and equa l to 5/3 in the Ke- 
plerian case (iPessah. Chan. & Psaltisll2006h . 

4.3. MRI Modes with Re > Rm 

Lets us consider the inviscid, poorly conducting limit de- 
scribed by = and r/ ^ 1. In this case, the marginally 
stable mode satisfying equation ( |70] i is given by 



and 



7n 



1 4 
?7 



(91) 



V^max value of the epicyclic 



kr — 



4-k2 



(87) 



The dependence of this critical wavenumber on the magnetic 
Reynolds number is shown on the left panel in Figure[3] which 
shows that for small magnetic Reynolds numbers kc oc Rm. 



As discussed in § 14.11 finding an analytic expression for 
the maximum growth rate and wavenumber associated with 
it is not as straightforward. The left panel of Figure |5] sug- 
gests that in the limit Re ^ oo and Rm ^ 1, the maxi- 
mum growth rate is linear in the magnetic Reynolds number, 
7inax oc Rm oc 77^^. This information can be used to derive 
asymptotic expressions for the dispersion relation jZSl ) and its 
derivative. The leading order contributions are given by 



In this case, 
frequency k. 

The dependence of both fcmax and 7inax in this hmiting case 
is shown with dashed lines in the left panels of Figures |4] and 
|5] respectively. The agreement between equations ( |90l ) and 
( I9TI 1 and the solutions to the full dispersion relation ( |25l l in 
the case i/ ^ and 77 ^ 1 is excellent, only breaking down 
close to magnetic Reynolds numbers of order unity. Note that 
even though the equations ( l90l ) and doTl l were derived under 
the assumption of an inviscid fluid, i.e., 1/ — 0, these expres- 
sions can describe the asymptotic behavior of both fc,„ax and 
7max for finite Reynolds numbers provided that the conditions 
Re ^ Rm and Rm ^ 1 are satisfied. 

Substituting the asymptotic expressions for /cmax and 7inax 
in equations ( |90] l and ( |9T| l into equation ( l54b we obtain the ra- 
tio between the amplitudes of the magnetic and velocity field 
perturbations 



bo 



(92) 



y2 

/ max 



+ 2k ?yA:^ax7max + K T] fc„,^^ + (k - 4)fc 



U2 

^max 



and 



-4 = 0, 



= 0, 
(88) 

(89) 



\/4 - k2 ' 

Therefore, inviscid, resistive MRl-unstable modes are domi- 
nated by magnetic field perturbations. Note that the ratio be- 
tween amplitudes increases linearly with resistivity. 

The asymptotic behavior for the angles characterizing ve- 
locity and magnetic field perturbations, equations ( fSOl l and 
(fSTl i. are given by 



respectively. 

Eliminating either jn 
and (|89]) we obtain 



or /cmax between equations 



and 



In the limit Re 



tan^v = 

tan6'b = 
00 and Rm 

lim e 



1 



2k^ 



V 



0, we obtain 

-0, 



4^2 



(90) 



lim ( 



'b = 77 



(93) 

(94) 

(95) 
(96) 
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Fig. 8. — Mean angle, Smcan = (6v + Sb) /2, defined by the fastest exponentia lly g rowi ng v elocity and magnetic perturbations for Keplerian rotation in various 
dissipative regimes. Note that Smcan = 7r/2 in the ideal MHD limit, see eqs. {83) and (84). Left: mean angle 9mcan as a function of the magnetic Reynolds 
number for different values of the Reynolds number The thick solid line denotes the inviscid limit, i.e.. Re oo. There is a critical magnetic Reynolds number 
Rriic < 1 that differentiates the asymptotic limits of Smoan for large and small Reynolds numbers. When the Reynolds number changes from Re <C Rm to 
Re 2> Rm, the mean angle evolves according to 9mean : [t^/'^ + arctan(K/2)]/2 — > 7r/2 for Rm > Rmc, while 9mcan : [■'r/2 + arctan(K/2)]/2 — > 7r/4 
for Rm < Rmc. Note that for Keplerian rotation [7r/2 + arctan(K/2)]/2 = 58° 17'. Middle: opening angle Sdiff as a function of the Reynolds number for 
different values of the magnetic Prandtl number. The thick solid line corresponds to Pm = 1. The thin solid lines to the right correspond to smaller values 
Pm = 10~^, 10~^, . . .. The thin solid lines to the left correspond to Pm = 10, lO'^, . . .. Note that for small Re, 9mcan ^ ■t/4 for Pm < 1, while 
fmcan — > [''"/2 + arctan{K/2)] /2 for Pm 2> 1. Right: mean angle (?mGan as a function of the magnetic Reynolds number for different values of the Reynolds 
number The thick solid line corresponds to the ideal conductor limit, i.e., Rm — > oo. The thin solid lines, in decreasing order, correspond to Re = 10, 1, . . .. 
For Reynolds numbers larger than unity, the mean angle is independent of Re regardless of the value of Rm. 



We therefore conclude that in the regime of large Reynolds 
numbers and small magnetic Reynolds numbers, magnetic 
field perturbations are larger than velocity field perturba- 
tions, both fields tend to be orthogonal and aligned with 
the azimuthal and radial directions, respectively, see Fig- 
uresEKH andlH 

4.4. MRI Modes with Re = Rm < 1 

When the magnetic Prandtl number is unity, and with rj = 
I' ^ 1, the marginally stable mode satisfying equation ( iTOl ) is 
given by 



vr- 



(97) 



The dependence of this critical wavenumber on the Reynolds 
number is shown in the middle panel in Figure|3] Incidentally, 
equation ( |97] i corresponds to the asymptotic limit 77 ^ 1 of 
equation (|87] |. 

It is not hard to see that the leading order contributions to 
the dispersion relation and its derivative in the limit Re = 
Rm ^ 1 are identical to the ones obtained in the case Rc 
00 and Rm ^ 1. Therefore , all the expressions for fcmax, 
7max, and dh, derived in § 14.31 are also valid in this case. 
The dependence of both fcmax and 7max in this limiting case 
is shown with dashed lines in the middle panels of Figures |4] 
and [ 5] res pectively. The agreement between equations (llOll i 
and ( IIO2I 1 and the solutions to the full dispersion relation jZSl l 
in the case h' ~ -q, i.e., Pm = 1, is also excellent in this case. 

4.5. MRI Modes with Re < Rm 

Lets us consider next the highly viscous, ideal conductor 
limit described by ^ 1 and r/ = 0. In this case, the 
marginally stable mode satisfying equation (fTOl i is given by 



kc — 



1/3 



(98) 



The right panel of Figure |5] suggests that in the limit Re ^ 
1 and Rm oo, the dependence of the maximum growth 



on the Reynolds number is 7max oc Re^^^ cx v^^/"^. This 
information can be used to derive asymptotic expressions for 
the dispersion relation dZSl ) and its derivative. The leading 
order contributions are given by 



,2^4 



)7iax + - 4)fci 



2 

max 



0, 



and 



max /max 



-4 = 0, 



(99) 



(100) 



respectively. 

Eliminating either 7n 
we obtain 



kvYisy- — 



between these equations 



(101) 



and 



O'max 



4- 



2VK 



(102) 



Eli minati ng the epicyclic frequency between equations ( I102l i 
and ( llOlb we obtain 

4 - v'^k'^ 



7^ 

/max 



^^^^max 



(103) 



The dependence of both fc,nax and 7,nax in this limiting case 
is shown with dashed lines in the right panels of Figures |4] 
and |5] respectively. The agreement between equations (llOll l 
and ( IIO2I 1 and the solutions to the full dispersion relation ( l25T l 
in the case v 1 and ?/ = is excellent, only breaking 
down close to Reynolds numbers of order unity. Note that 
even though the equations ( llOlb and ( I102l i were derived un- 
der the assumption of a perfectly conducting fluid, i.e., rj = 0, 
these expressions can describe the asymptotic behavior of 
both fcmax and 7max for finite Reynolds numbers provided that 
the conditions Re <C Rm and Re ^ 1 are satisfied. 

Substituting the asymptotic expressions for /cmax and 7inax 
in equations (llOlb and ( 1102b into equation (|54] | we obtain the 
relative amplitude of the magnetic and velocity field perturba- 
tions 

= ^. (104) 
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Fig. 9. — Geometrical representation of viscou.s. resistive MRI modes for varying Reynolds and magnetic Reynolds numbers. The black/gray lines denote 
velocity/magnetic components of the most unstable mode. 



Therefore, viscous, conducting MRI-unstable modes are also 
dominated by magnetic field perturbations. In this case, the 
ratio between amplitudes increases only with the square root 
of the viscosity. 

The asymptotic behavior for the angles characterizing ve- 
locity and magnetic field perturbations, equations ( l50b and 
( fSTT ). are given by 



tan6'v = — 



and 



In the limit Re 



tanPh = 
and Rm 
lim 6v 

V — >OQ 



V'2r/K(4 - k2) 



i/(4 - K^) 
2k 

> oo, we obtain 



: arctan 



lim 9\, — — . 



(105) 

(106) 

(107) 
(108) 



For Keplerian rotation, the angle between the fastest grow- 
ing velocity field perturbation and the radial direction is given 
by 9^ = arctan(l/2) = 26°34'. Therefore, the opening an- 
gle between the planes containing velocity and magnetic field 
perturbations is 6'diff ^ Oy, — 9^ = 63°26' and their mean 



value is Oy, 



7diff 

(0b 



+ 9^) /2 ^ 58° 17'. The right panels of 



Figures I?] and [8] show that both of these results are in perfect 
agreement with the asymptotic behavior of the full solutions 
derived directly from the original dispersion relation ( |25] l. 

We therefore conclude that in the regime of small Reynolds 
numbers and large magnetic Reynolds numbers, magnetic 
perturbations are larger than velocity perturbations. In this 
case, however, the perturbed magnetic and velocity fields are 
not orthogonal. The perturbed magnetic field tends to be 
aligned with the azimuthal direction but the velocity field per- 
turbations do not tend to be aligned with the radial direction. 
The angle between both fields is determined entirely by the 
epicyclic frequency k, see Figures |7]|8] and|9] 

4.6. MRI Modes with Re = Rm > 1 

As we show in Appendix lAl when the Reynolds and mag- 
netic Reynolds numbers are large enough the solutions to the 
dispersion relation ( |25] | tend smoothly to the solutions of the 
dispersion relation (l72l i for the idealized case. Therefore, in 
this limit we recover all the expressions derived in 14.21 

5. PHYSICS OF MAXIMALLY UNSTABLE MRI MODES 

We have shown that the expressions derived for the most 
unstable wavenumber, fcmax, and its associated growth rate, 
7max, obtained from the simplified equations dSSb and ( |99] l are 
good approximations to the solutions obtained directly from 
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the dispersion relation dZSl ). in the limits Re 3> Rm and Re <C 
Rm, respectively. We can now identify the various terms in 
the original set of equations of motion (iTt-dTOli that lead to 
equations dSST l and ( |99] l. This allows us to better understand 
the forces that act to destabilize magnetized fluid elements. 

For the sake of clarity we write the equations in this section 
with physical dimensions. We represent the temporal deriva- 
tives with a dot and the derivatives with respect to the vertical 
coordinate z with a prime. 

5.1. MRI Modes with Re > Rm 

The equations of motion that lead to the dispersion relation 
are given by 



5Vr = 2rin 



SBL 



= -{2-q)noSvr + -^SB' 

4:TTp ^ 



3" 



SBr = BMr + vSB'; 

SB^^-qrioSBr + ijSB';. 



(109) 
(110) 

(111) 

(112) 



Therefore, maximally unstable modes with Re ^ Rm are 
characterized by motions with radial acceleration due to the 
Coriolis force acting on azimuthally displaced fluid elements 
and magnetic tension. Azimuthal force balance is attained 
via the joint action of the Coriolis force acting on radially 
displaced fluid elements, radial advection of background flow, 
and magnetic tension in the azimuthal direction. The rate of 
change of the radial magnetic field perturbations is due to the 
creation of radial field out of the vertical background frozen 
into the radial velocity field with a vertical gradient and field 
diffusion. The rate of change of the azimuthal magnetic field 
perturbations is due to the shearing of radial magnetic field 
perturbations and field diffusion. 

5.2. MRI Modes with Re < Rm 



The set of equations that lead to the dispersion relation 
are given by 



(113) 



B, 



= -{2-q)noSvr + ^6B'^ + i^6v';, (114) 

(115) 
(116) 



SBr^BJvl , 
SBcf, = —qfloSBr 



In this case, maximally unstable modes with Re ^ Rm 
are characterized by fluid displacements that take place un- 
der force balance in both the radial and azimuthal directions. 
The Coriolis force acting on azimuthally displaced fluid ele- 
ments is balanced by the viscous force in the radial direction. 
Azimuthal force balance is attained via the joint action of the 
Coriolis force acting on radially displaced fluid elements, ra- 
dial advection of background flow, magnetic tension, and the 
viscous force in the azimuthal direction. The rate of change of 
the radial magnetic field perturbations is due to the creation of 
radial field out of the vertical background frozen into the ra- 
dial velocity field with a vertical gradient. Finally, the rate of 
change of the azimuthal magnetic field perturbations is due to 
the shearing of radial magnetic field perturbations. 



6. NON-IDEAL MRI-DRIVEN STRESSES AND ENERGY DENSITIES 

The i, j components of the mean Reynolds and Maxwell 
stresses associated with the velocity and magnetic field per- 
turbations are given by 

R,j{t) = {6v,{z,t)Svj{z,t)), (117) 
M,j{t) = {Sh{z,t)Sb,{z,t)) , (118) 

where the brackets denote mean values obtained via integra- 
tion over the disk scale-height, 2H. These mean values can 
also be calculated directly from the pert urbations in Fourier 
space according to (iPessah. Chan. & Psa ltis 2006) 



R,,it) = 2j2 ne[6v,{kn,t)6i*ikn,t)], (119) 

n=l 

oo 

M„(t) = 2^ TZe[Sk{k,„t)6b*{kn,t)]. (120) 

n=l 

where Tle[] stands for the real part of the quantity be- 
tween square brackets. The off-diagonal components of these 
stresses relate to angular momentum transport while their 
traces relate to the kinetic and magnetic energy in the per- 
turbations. 

At late times, during the exponential growth of the instabil- 
ity, the branch of unstable modes will dominate the growth of 
the perturbations and the most important (secular) contribu- 
tion to the mean stresses will be given by the most unstable 
mode. The leading order contribution to these stress com- 
ponents are thus obtained by considering the most unstable 
solutions to the set of equations (TTll-lfTOll. which are given by 
equation ( l63T l when k = fcmax with 7+ 



7n 



6.1. Non-ideal MRI Stresses 

A measure of the angular momentum transport driven by 
the most unstable MRI modes and mediated by the correlated 
perturbations in the velocity and magne tic fie lds is obtained 
by setting (i, j) = (r, 0) in equations (II 17b and jl 18b . To 
leading order in time we obtain 



1 v^e^-'" 



2 vl+bl 



1 62e27„ 



2 «2+62 



sin(26'v 



sin(26'b) . 



(121) 
(122) 



The results derived in §21 together with Figures |2l[8j and|9l 
show that the angles 9^ and 9\, corresponding to the most un- 
stable mode, fcmax, always satisfy 



0<6'v < arctan (^-j < - , 

-<9h< — . 
2-4 



(123) 
(124) 



Both of these inequalities show explicitly that the mean 
Reynolds and Maxwell stresses will be, respectively, positive 
and negative. 



Rr4>it) > 



and 



Mr4,it) < 0. 



(125) 



This, in turn, implies that the mean total MRI-driven stress 
will be always positive, i.e.. 



Tr4,it) ^ Rr4,{t) - Mr4,{t) >0, 



(126) 
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Fig. 10. — Ratio between the Maxwell and the Reynolds stresses, —M^^/R^^, for Keplerian rotation in different dissipative regimes. Left: Ratio between 
the Maxwell and the Reynolds stresses as a function of the magnetic Reynolds number for different values of the Reynolds number. The thick solid Hne denotes 
the inviscid limit, i.e.. Re oo. The thin solid Hnes, in increasing order according to their asymptotic value at high magnetic Reynolds numbers correspond to 
Re = 1, 0.1, . . .. For Rm > 1, the ratio is independent of Rm regardless of the value of Rc. Middle: Ratio between the Maxwell and the Reynolds stresses as a 
function of the Reynolds number for different values of the magnetic Prandtl number. The various curves, from left to right, correspond to Pm = 10'^, 10'^, .... 
The thick solid line corresponds to the case Pm = 1. Right: Ratio between the Maxwell and the Reynolds stresses as a function of the Reynolds number 
for different values of the magnetic Reynolds number. The thick solid line con'esponds to the ideal conductor Hmit, i.e., Rm oo. The thin solid lines, in 
increasing order according to their asymptotic value at high Reynolds numbers correspond to Rm = 1, 0.1, 0.01. The Maxwell stress is larger than the Reynolds 
stress for any combination of the Reynolds and magnetic Reynolds numbers. The minimum value of this ratio is achieved in the ideal MHD regime, where 
— M^^/ R^^ = 5/3, for Kepleiian rotation. 



driving a net outward flux of angular momentum for any com- 
bination of Reynolds and magnetic Reynolds numbers. 

We conclude this section by calculating the ratio 
—Mrc/,{t)/Rr^{t) at late times during the exponential growth 
of the instability. We obtain 



lim 



-Mr4>{t) _ vl sin(26'v) 



bl sin(20b) 



(127) 



Using the definitions for the ratio wo/feo (eq- ll54l ) and the 
angles 9^ and 0\y (eqs. llSOl and ||5TI \ it can be seen that the 
magnitude of the Maxwell stress, ~AIr^{t), is always larger 
than the magnitude of the Reynolds stress, Rr^{t), provided 
that the flow is Rayleigh-stable, i.e.. 



Mr4,{t) > Rr^{t) for < g < 2 . 



(128) 



Figure [TO] shows the ratio between the Maxwell and the 
Reynolds stresses in various dissipative regimes for Keple- 
rian rotation. Note that when the Reynolds and magnetic 
Reynolds numbers are large enough we recover the result 
—Mr4,/Rr4, — 5/3, which coincides with the value of th is 
ratio in the ideal MHD case (iPessah. Chan. & Psaltisll2006h . 

6.2. Non-ideal MRI Energetics 

The mean energy densities associated with the perturba- 
tions in the velocity and magnetic field are given by 



[Rr 



(129) 
(130) 



Substituting the expressions for the most unstable MRI-driven 
perturbations from equation (|63] | into the definitions for the 
diagonal components of the Reynolds and Maxwell stresses, 
equations ( II 17l )-( ll 18b . respectively, we obtain 



EK{t)-- 



EM(t)-- 



1 vy^" 



2 + 

1 62e2^max« 



(131) 



(132) 



Using the definitions for the ratio vo/bo, equation ( |54l ). it 
can be seen that for non-ideal MRI modes the mean energy 
associated with magnetic perturbations is always larger than 
the mean energy corresponding to velocity perturbations as 
long as the flow is Rayleigh-stable. Figure [TT] shows the ra- 
tio between the magnetic and the kinetic energy densities in 
various dissipative regimes for Keplerian rotation. 

In the limit of large Reynolds and magnetic Reynolds num- 
bers, we recover the result Em /Ek = 5/3, which coincides 
with the value of this ratio in the ideal MHD case and also 
with the ratio between the magnitu des of the Maxwell and 
Reynolds stresses in the ideal case (iPessah. Chan. & PsaltisI 
|2006) . Note, however, that for arbitrary Reynolds and mag- 
netic Reynolds numbers, it is no longer true that the ratio be- 
tween mean magnetic and mean kinetic energies is equal to 
the ratio between the magnitude of the mean Maxwell and 
the mean Reynolds stresses. This can be seen by comparing 
Figures[IO]and[n] 

Finally, comparing equations (11211 ) and ( I122l i with ( I131l l 
and ( II32I 1. it immediately follows that 



Rr4>{t) < Exit) , 
-Mr4t)<EM{t). 



(133) 
(134) 



This result, in turn, implies that the total mean energy associ- 
ated with the perturbations, E{t) — EK{t) + EM{t), sets an 
upper bound on the total mean stress, i.e.. 



for any Reynolds and magnetic Reynolds numbers. 



(135) 



7. SUMMARY & DISCUSSION 

We investigated the effects of viscosity and resistivity on 
the stability of differentially rotating plasmas threaded by a 
magnetic field perpendicular to the shear. We have shown 
that the most powerful incompressible MRI modes are exact 
solutions of the MHD equations for arbitrary combinations of 
the Reynolds and magnetic Reynolds numbers. We have de- 
rived analytical expressions for the eigenfrequencies as well 
as for the eigenmodes describing the MRI in viscous, resis- 
tive media and provided a detailed description of the physical 
properties of these modes in various dissipative regimes. 
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Fig. 1 1 . — Ratio between the magnetic and the kinetic energy densities contained in MRI-driven perturbations, Em / Ek, for Keplerian rotation in different 
dissipative regimes. Left: Ratio between the magnetic and the kinetic energy densities as a function of the magnetic Reynolds number for different values of the 
Reynolds number The thick solid line denotes the inviscid limit, i.e., Rc — > oo. The thin solid lines, in increasing order according to their asymptotic value at 
high magnetic Reynolds numbers correspond to Re = 1, 0.1, . . .. For magnetic Reynolds numbers larger than unity, the ratio is independent of Rm regardless 
of the value of Rc. Moreover, the asymptotic value of this ratio for Rm <g; 1 is independent of the Reynolds number. Middle: Ratio between the magnetic 
and the kinetic energy densities as a function of the Reynolds number for different values of the magnetic Prandtl number The various curves, from left to 
right, correspond to Pm = 10^, 10, . . .. The thick solid line corresponds to the case Pm = 1. Right: Ratio between the Maxwell and the Reynolds stresses 
as a function of the Reynolds number for different values of the magnetic Reynolds number The thick solid line coiresponds to the ideal conductor limit, i.e., 
Rm — > oo. The thin solid lines, in increasing order according to their asymptotic value at high Reynolds numbers correspond to Rm = 1,0.1, .. .. Unlike 
the ratio between stresses, the ratio between energy densities seems to decrease monotonically with Rc for any value of the magnetic Reynolds number The 
magnetic energy density is larger than the kinetic energy density for any combination of the Reynolds and magnetic Reynolds numbers. The minimum value of 
this ratio is achieved in the ideal MHD regime, where Ej^j /Ex = 5/3, for Keplerian rotation. 



We have shown that the scalings derived for the marginally 
stable mode, the most unstable wavenumber, and the 
maximum growth rate with magnetic Reynolds number, 
fcc, fcmax, 7max oc Rm, vaHd for resistive, i nviscid plasmas, 
see eq uations dSTl i. ( l90l l and doTl i. as well as Sano & Miyamal 
1 19991 also hold when finite Reynolds numbers are involved. 
This is true as long as the magnetic Prandtl number is of order 
unity or smaller, as it is usually the case in many astrophysical 
systems (such as accretion disks around cataclysmic variables 
and young stellar objects, as well as the Sun) and also in MRJ 
laboratory experiments. Furthermore, we have addressed in 
detail, for the first time to our knowledge, the physical prop- 
erties of the MRl in highly viscous, slightly resistive media. 
These conditions are expected to be found in the hot, diffuse 
gas in galaxies and galaxy clusters. In this case, we found 
that the critical wavenumber for the onset of the MRl, the 
most unstable wavenumber, and the maximum growth rate 
scale with the Reynolds and magnetic Reynolds numbers ac- 
cording to fee OC (ReRm)^/'^ and fcinax,7max oc Re^/^, see 



equations i 



fToD ) and (fToa . 



We have provided a thorough geometrical description of the 
viscous, resistive MRl modes in terms of the angles that define 
the planes containing the velocity and magnetic field pertur- 
bations. In the ideal MHD limit, these planes are orthogonal, 
with the plane containing the velocity disturbances laying at 
45° with respect to the radial direction. We have shown that 
velocity and magnetic field perturbations are still orthogonal 
if the magnetic Prandtl number is unity, but that the planes 
containing them tend to be aligned with the radial and az- 
imuthal directions, respectively, when the Reynolds number 
increases. In the regime of large Reynolds numbers and small 
magnetic Reynolds numbers, magnetic and velocity field per- 
turbations tend to be orthogonal and aligned with the az- 
imuthal and radial directions, respectively. On the other hand, 
in the regime of small Reynolds numbers and large magnetic 
Reynolds numbers, the perturbed magnetic field tends to be 
aligned with the azimuthal direction but the velocity field per- 
turbations do not tend to be aligned with the radial direction. 
The angle between both fields is determined entirely by the 
epicycUc frequency k. It would be very interesting to un- 



derstand to what extent this geometrical dependence of MRl 
modes on the Reynolds and magnetic Reynolds numbers in- 
fluences the physical properties of kinetic and magnetic cells 
in fully developed viscous, resistive MHD turbulence. 

In the ideal MHD limit, the exact (primary) MRl modes 
are known to be unstable to parasitic (secondary) instabili- 
ties dGoodman & Xu|[l994l) . These parasitic modes have long 
been suspected to enable the mechanisms that disrupt the pri- 
mary modes providing an avenue toward saturation. The so- 
lutions derived in this work describe the dynamics of primary 
MRl modes in viscous, resistive media enabling the study of 
parasitic instabilities for arbitrary combinations of Reynolds 
and magnetic Reynolds numbers. The modifications in the 
relative directions of the velocity and magnetic field perturba- 
tions characterizing the primary viscous, resistive MRJ modes 
described above can have an important impact on the devel- 
opment and evolution of parasitic instabilities in the presence 
of dissipation. 

We have shown that, for any combination of the Reynolds 
and magnetic Reynolds numbers, the mean Reynolds stress, 
Rr4> — {Svr{z, t) 5v^{z, t)) , is always positive and the mean 
Maxwell stress, M^^ — {Sbr{z, t) 5b^{z, t)), is always nega- 
tive. This implies that the mean total stress, Trcf, = Rr4, — Mr4, 
is always positive, leading always to an outward transport of 
angular momentum. We have also demonstrated _that both 
the ratio between magnetic and kinetic stresses, — Mr^/i^r^, 
and the ratio between magnetic and kinetic energy densities, 
Em I Ek, are always dominated by the magnetic contribution. 
These last two statements, support a somewhat unexpected 
result since it is tempting to think that velocity perturbations 
would dominate both the transport of angular momentum and 
the energy density in highly resistive, inviscid plasmas. It 
would be very interesting to understand if and how the value 
of these ratios in the saturated turbulent state vary as a func- 
tion of the Reynolds and magnetic Reynolds numbers. 
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Fig. 12. — Mean total stress responsible for angular momentum transport, = R^^ — Mj.^, calculated according to equations U2U and (122), in various 
dissipative regimes for Keplerian rotation. Left: Mean total stress TV ^ as a function of the magnetic Reynolds number for different values of the Reynolds number 
The thick solid line denotes the inviscid limit, i.e.. Re oo. The thin solid lines, in decreasing order, correspond to Rc = 10, 1, . . . , 10~^. For magnetic 
Reynolds numbers larger than unity, this ratio is independent of Rm regardless of the value of Rc. Moreover, the asymptotic value of this ratio for Rm 1 is 
independent of the Reynolds number Middle: Mean total stress f^^ as a function of the magnetic Prandtl number for different values of the Reynolds number 
From left to right, the curves con'espond to Pm = 10'^, 10^, . . . , 1 (thick solid line), . . . , 10"''. Right: Mean total stress T^^ as a function of the Reynolds 
number for different values of the magnetic Reynolds number The thick solid line corresponds to the ideal conductor limit, i.e., Rm oo. The thin solid lines, 
in decreasing order, correspond to Rm = 10, 1, 0.1. Note that the stress scale in these plots is arbitrary. 



for inv iscid, resistive MHD. The simulations in lSano & Stond 
(l2003h show that for small magnetic Reynolds numbers the 
stresses at saturation increase rapidly with increasing Rm 
and that there exists a critical magnetic Reynolds number, 
of order unity, beyond which turbulent stresses are rather 
insensitive to Rm. The fact that this same behavior is indeed 
seen when the stresses are due to viscous, resistive MRl 
modes (see Figure [T2] | rises the question of how strong is 
the influence of long-lived, channel-like modes on the fully 
developed turbulent state reached in shearing box simulations 
with net magnetic flux through the vertical boundaries. 

Systematic numerical studies of viscous, resistive MHD 
shearing flows have begun to uncover the dependencies of 
microphysical dissipation on the mean transport properties 
of MRJ-driven turbulence. Numerical simulations with both 
zero fer omang et al. 2007) and non-zero net magnetic fluxes 
(iLesur & Longarettiil2007l) lead to the conclusion that angu- 
lar momentum transport increases with increasing magnetic 
Prandtl number when the Reynolds number is held constant. 
This behavior can also be identified when examining the 
stresses due to viscous, resistive MRl modes, see Figure [12] 

The effects of varying the Reynolds number at fixed mag- 
netic Prandtl number on the saturation of MRI-driven turbu- 
lence are currently rather uncertain (see, in particular, the dis- 
cussion in Lesur & Longaretti i2007b . The global trends ex- 
hibited by the available simulations suggest that the stresses at 
saturation increase with increasing Reynolds number for mag- 
netic Prandtl numbers smaller than unity while the opposite 
behavior is observed for magnetic Prandtl numbers larger than 
unity. If confirmed, these results suggest that the mechanisms 
leading to saturation might operate differently depending on 
whether the magnetic Prandtl number is larger or smaller than 
unity. In any case, having obtained a better understanding 
of the behavior of the most unstable MRl modes as a func- 



tion the magnetic Prandtl number it would be very interest- 
ing to follow the evolution of the viscous, resistive MRl from 
the linear to the non-linear regime. By performing numer- 
ical simulations with the same Prandtl number (both larger 
and smaller than unity) and different Reynolds numbers we 
could see whether there is an inversion of the trends observed 
in the linear regime (i.e., higher stresses at higher Reynolds 
numbers for fixed magnetic Prandtl numbers) after the exact 
solutions break down. The comparison between the late time 
behavior of the viscous, resistive MRl modes and fully devel- 
oped MHD turbulence with dissipation will shed light into the 
mechanisms that lead to the saturation of the MRl in different 
dissipative regimes. 

Finally, most current numerical algorithms employ finite 
difference methods (with constrained transport for the evo- 
lution of the magnetic field). The leading order errors in 
first-order upwind methods behave like diffusion, however, 
in second-order central difference methods, the leading or- 
der errors are dispersive (artificial viscosities are usually em- 
ployed to damp unphysical oscillations near shocks). The 
comparison between numerical solutions from different al- 
gorithms can provide estimates of these errors. However, it 
is difficult to quantify these numerical artifacts based on an- 
alytical studies of ideal MHD. The analytical solutions de- 
rived in this paper, on the other hand, describe the effects of 
arbitrary combinat i ons o f viscosity and resistivity [see also 
iLesaffre & BalbusI (l2007h who derived results to leading or- 
der in [r] — v)/k'\. It should now be possible to better mea- 
sure the numerical viscosity and resistivity for a wide range of 
Reynolds and magnetic Reynolds numbers by comparing nu- 
merical solutions and analytical solutions of non-ideal MRl. 
The results presented in this paper provide ideal benchmarks 
to the study numerical artifacts generated by different algo- 
rithms in various dissipative regimes. 



APPENDIX 



ANALYTICAL SOLUTION TO QUARTIC EQUATION AND IDEAL MHD LIMIT 



The solutions to a depressed quartic equation of the form 

+ aa'^ + 13a + \^ 0. 



(Al) 
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are given by^ 



= ±a^ 



where y is any of the solutions of the cubic equation 



3a y 13/ A 

— + - ±6 



" , y 



a3 _ oA _ 
"2 2 8 



0. 



(A2) 



(A3) 



In the special case fi 0, the solutions to equation ( IA3I) take simple forms. In order to see how the general solutions reduce 
to the simple cases, we write the cubic equation as 



y+2) [(2^ + ") 



It is easy to see that, iiy ^ —a/2. 



(3/4 



V«/4 + y/2 

Using the above identity, the general solution iA2i becomes 



A/(y + a)2-A. 



/3/4 



A 



where we have defined 



If we choose the root so that 



it is easier to take the limit 



A 



So: y 
T ^ 2 



and A = \/{y + a)'^ - X . 



lim V = . 

/3^0 ^ 2 



where Aq and Aq are defined in equations (TmI i and ( |75] l. 



lira (7 lim^ ± V-A ± VA = ±\/ -| T Y ^ - A = ± V -Ao ± = ctq 



(A4) 
(A5) 

(A6) 
(A7) 
(A8) 

(A9) 
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